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Abstract 

We consider the propagation of sine-Gordon kinks in a planar curved 
strip as a model of nonlinear wave propagation in curved wave guides. The 
homogeneous Neumann transverse boundary conditions, in the curvilinear 
coordinates, allow to assume a homogeneous kink solution. Using a simple 
collective variable approach based on the kink coordinate, we show that 
curved regions act as potential barriers for the wave and determine the 
threshold velocity for the kink to cross. The analysis is confirmed by 
numerical solution of the 2D sine-Gordon equation. 

PACS numbers: 74.50. +r, 05.45.Yv, 85.25. Cp 

Recent advances in micro-structuring and nano-structuring technology have 
made it possible to fabricate various low-dimensional systems with complicated 
geometry. Examples are photonic crystals with embedded defect structures 
such as microcavities, wave guides and wave guide bends p^; narrow construc- 
tions (quantum dots and channels) formed at semiconductor heterostructurcs 
0, magnetic nanodisks, dots and rings e ^ c - 

It is well known that the wave equation subject to Dirichlet boundary con- 
ditions has bound states in straight channels of variable width [Sj and in curved 
channels of constant cross-section [5]. Spectral and transport characteristics 
of quantum electron channels [7] and wave guides in photonic crystal [Sj are 
essentially modified by the existence of segments with finite curvature. 

Until recently there have been a few theoretical and numerical studies of the 
effect of curvature on properties of nonlinear excitations. The dynamics of a ring 
shaped Josephson fluxons and their collisions was studied in 111)1 ITT] . Non- 
linear whispering gallery modes for a nonlinear Maxwell equation in microdisks 
were investigated in , the excitation of whispering-gallery-type electromag- 
netic modes by a moving fluxon in an annular Josephson junction was found 
in 13 . Nonlinear localized modes in two-dimensional photonic crystal wave 
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guides were studied in ^3] . A curved chain of nonlinear oscillators was consid- 
ered in |15| and it was shown that the interplay of curvature and nonlinearity 
leads to a symmetry breaking when an asymmetric stationary state becomes 
energetically more favorable than a symmetric stationary state. Propagation 
of Bose-Einstein condensates in magnetic wave guides was found quite recently 
in Single- mode propagation was observed along homogeneous segments of 
the wave guide while geometric deformations of the microfabricated wires lead 
to strong transverse excitations. 

The aim of this article is to study the motion of sine-Gordon (sG) solitons 
moving in a two dimensional finite domain. Specifically we treat a planar curved 
wave guide whose width is much smaller than its entire length. We consider 
homogeneous Neumann boundary conditions (zero normal derivative) on the 
boundaries of the domain. Using a simple collective variable analysis based 
on the kink position we show that a region of non zero curvature in a wave 
guide induces a potential barrier for the wave as it can be verified in the 2D 
simulations made by Femlab finite element software. This is different from the 
case of transverse Dirichlet boundary conditions where studies on the (linear) 
Schrodinger equation show the existence of a localized mode which will trap 
waves in the curved region . 

In the following we study the sine Gordon equation as a model for nonlinear 
wave propagation in planar curved wave guides. A physically relevant example 
is a Josephson junction constructed of two straight segments joined by a bent 
section (see Fig. 1). In two spatial dimensions and disregarding the effects of 
loss and external power inputs the sG equation reads 
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together with boundary conditions and two initial conditions. For a Josephson 
junction the function (f>{x, y, t) is the normalized phase difference of the Cooper 
pair wave functions across the insulating barrier. 
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Figure 1: The bend space for the 2+1 D sG equation. 

To solve equation Q in the special domain shown in Fig. ^ we transform 
the spatial coordinates [x, y) into (s, u), where s is arc length of the central line 
and u is the transversal coordinate orthogonal to s. A simple geometry with 
a curved section, D, of constant radius R has been considered here, where the 
curvature is given by 

7 (») = ( ~R> When *i< s < s 2- (2 ) 

I 0, elsewhere 
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We parameterize the centerline according to (x,y) — a(s) — [a(s),b(s)] 
and choose the normalized parameterization such that ||cr'(s)|| = 1. The unit 
normal vector n then becomes [-b'(s), a' (s)]. Considering for simplicity wave 
guides with the width smaller than the radius of the curvature R > w/2, we can 
describe the wave guide domain by means of the following set of coordinates 
(x, y) = [a(s) - ub'(s),b(s) + ua'(s)] as in Ref. 0. 

By i we denote the total length of the strip and < s < £. The curved 
region D corresponds to the interval si < s < S2- The width of the strip is w so 
that —w/2 < u < w/2. We can determine a function 8(s) such that the tangent 
vector in the curved region can be written in the form t={cos[6>(s)], sin[6*(s)]}, 
where 8(s) = (n/2)(s — si)/(s2 — si). The signed curvature 7(s) now becomes 
7(s) = a'(s)b"(s) - a"{s)b'(s) = 9'(s). The solution space [0;£] x [-w/2; w/2] 
we denote by Q. It is Ricmannian with the metric gij, where i and j are s 
or u. The determinant of the Jacobian matrix of the variable change leads to 
.9 = 1- uj(s). 

The Lagrange function of the sine Gordon equation in Cartesian coordinates 
(x,y) can be transformed into the new variables (s,it) 



1 J.2 ^ J.2 1 1 2 



2 Yt 2g 2^s 2 Vu 

[1 — cos(0)] \gduds . (3) 



Variation of the above Lagrange function provides the sG equation for <fi in the 
coordinate space (s, u, t), 

d 2 cj) 1 d fld(j>\ 1 d ( d<j>\ . n 

To this we add the Neumann boundary conditions corresponding to a Josephson 
junction with no magnetic field or external current 

— = 0, on the boundary dil. (5) 

Both cqs. Q and 10} can be solved numerically in their respective domains. 
However first we investigate kink propagation for the sG equation (@J) in the 
curvilinear coordinate space using a simple collective coordinate approach for 
the kink position. We assume that our solution <p only depends on s and t, 
neglecting the variation along the it-direction, 4> u = 0. We denote this solution 
as <fio=(f>o(s,t) and it reads 

4>o(s,i) = 4arctan e s " s(t) , (6) 

where S(t) is the position of the center part of the wave. 

The evolution of S is given by the reduced Lagrangian obtained by substi- 
tuting (f>Q in eq. JBJl and performing the spatial integrations over s and u. We 
have L = L{(/)q}. The reduced Lagrangian is then 

dS K 2 



L = 4w[ — ) -U(S), (7) 
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where the effective potential energy U(S) is 

1 



U 



^4>l + [1 - cos(0)] 



gduds 
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1 - w/(2R) J S1 
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sinh(s 2 - si) 



6R 2 cosh(5' — si) cosh(s2 — S) 



(8) 



From this expression it is clear that the term U and therefore the curved section 
acts as a potential barrier for the kinks. The maximum of the barrier is w 3 /3i? 2 . 
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Figure 2: The wave is transmited for R = 2, w — 2 and = 0.25, while 
v c w 0.29. 

The initial kinetic energy of the kink, which is located in the straight section, 
far away from the curved region and from the domain edge to avoid interference 
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from the boundaries, with initial velocity vq, is 

r w/2 ^ 



w/2 



<j/ t du ds = 4wv . 



(9) 



When To is large enough to overcome the potential U the wave crosses the barrier 
running along the channel as shown in the contour-plot sequence of Fig. [21 The 
minimal velocity for the kink to cross the barrier is calculated by comparing the 
kinetic energy Tq to the potential energy U . It is approximately 



12R 



tanh 



S2 - Si 



(10) 



On the other hand when the wave is not initially fast enough, the potential 
barrier reflects the wave and it returns with a negative velocity (see Fig. EJ. 
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Figure 3: The wave is reflected for R = 2, w = 2 and i>o = 0.25, while v c s» 0.29. 

We now proceed to compare quantitatively the predictions of this simple 
theory with the solution of the 2D problem. We do this by plotting the phase 
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space (S, dS/dt) both for the collective variable dynamics (energy levels) and 
the kink position S and velocity dS/dt estimated from the 2D solution. The 
initial position of the kink is fixed well away from the curved region 20 < 
S < 20 + Rtt/2. Then the initial velocity of the soliton determines completely 
its trajectory as can be seen in the phase space (S,dS/dt) shown in Fig. 0] 
When the initial velocity is lower than vq w 0.2 the orbits do not cross the 
curved region. On the other hand for higher initial velocities the orbits span the 
complete circuit through the S axis. Due to the existence of effective potential 
barrier the velocity of kink decreases in the curved segment. Fig. 01 shows that 
the orbits calculated from the numerical solution of the 2D system are are in 
very good agreement with the ones given by the collective coordinate approach. 
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Figure 4: Phase space (S,dS/dt) of the kink position and velocity for R = 2 
and w — 2. Comparison of simulations (continuous lines) and the collective 
coordinates (dashed lines). 

The simple analytical approach developed above provides a good description 
of the kink dynamics for small velocities, otherwise the ansatz used in Eq. (JJjJ 
may be too crude for large ones. In that case the kink ansatz should also 
depend on the variable u as the contour plot of Fig. [S] shows. The distance in 
s covered by the wave increases as u decreases from w/2 to —w/2. Thus the 
soliton velocity depends on the level at u axis and Lorenz contraction induces 
the different kink widths observed for different values of u. 




Figure 5: Contour plot of the wave at the curved section for R = 2, w = 2 and 
v = 0.20. 
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The numerical simulations are done in the cartesian coordinates using the 
Femlab finite element software to define the domain and discretizing it by 
a triangular mesh. The time evolution has been performed with a differential 
equation solver. 

We consider a nonlinear wave propagation in a curved planar wave guide 
using as a model kink solutions for the sine-Gordon equation. Specifically the 
domain considered f2 in the plane is made of two rectangular regions joined 
by a bent section of constant curvature. The transverse homogeneous Neuman 
boundary conditions allow us to consider a homogeneous kink solution. Follow- 
ing this we develop a simple collective variable theory based on the kink position. 
This shows that curved regions act as potential barriers for the waves in contrary 
to the case of Dirichlet boundary conditions. We calculate the treshold velocity 
for the kink to cross and it is in excellent agreement with the solutions of the 
full problem calculated with Femlab finite clement program. The phase-space 
of the system matches also well the one obtained from the numerical solution. 

This study shows that it is possible to trap fluxons in a Josephson junction, 
in the region between two successive bends. Choosing a conveniently small 
damping term, one may decrease the kinetic energy of the soliton as it comes 
into this region so that it will stay there. This feature could be applied to 
electronic devices for storing binary data. 
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